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Abstract 

We investigate the nonuniform motion of a straight screw dislocation in infinite 
media in the framework of the translational gauge theory of dislocations. The equa- 
tions of motion are derived for an arbitrarily moving screw dislocation. The fields 
of the elastic velocity, elastic distortion, dislocation density and dislocation current 
surrounding the arbitrarily moving screw dislocation are derived explicitely in the 
form of integral representations. We calculate the radiation fields and the fields 
depending on the dislocation velocities. 
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1 Introduction 

For many years, the investigation of the nonuniform motion of dislocations has attracted 
the interest of researchers in different fields such as physics, material science, continuum 
mechanics, seismology and earthquake engineering. Usually, the motion of dislocations 
is investigated in the framework of incompatible elastodynamics where the dislocation 
density tensor and the dislocation current tensor are given as source terms of the elastic 
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fields (see, e.g., Mura [1, 2], Lardner [3], Giinther [4]). The behaviour of the motion of 
a dislocation is somehow particular, because at any time the fields are determined not 
only by the instantaneous values of the velocity (or higher derivatives of the position with 
respect to time), but also by the values of these quantities in the past [5, 6]. As Eshelby 
[5] succinctly put it: 'The dislocation is haunted by its past'. This fact is based on the 
property that Huygens' principle is not valid in two dimensions (see, e.g., [7]). Solutions 
for nonuniformly moving screw dislocations have been given by Nabarro [8], Eshelby 
[5, 6], Kiusalaas and Mura [9], and Markenscoff [10]. Results for the nonuniform motion 
of a ghding edge dislocation have been given by Kiusalaas and Mura [11], Markenscoff 
and Clifton [12], and Brock [13]. 

Immediately, the question comes up, what is the suitable theory for an improved 
approach of a continuum theory of dislocation dynamics. Generalized theories of elasticity 
are, for instance, nonlocal elasticity [14, 15], gradient elasticity [16, 17] and dislocation 
gauge theory [18, 19, 20, 21]. However, in the dynamical version of nonlocal elasticity 
it is not obvious which dynamical kernel should be used. In the dynamical extension of 
gradient elasticity [22, 23] the choice of 'dynamical' gradients and length scales is based 
on ad-hoc assumptions. On the other hand, dislocation gauge theory is not based on 
such assumptions. The choice of the dynamical state variables is given by a canonical 
field theoretical framework [21]. For this reason, dislocation gauge theory seems to be the 
appropriate theory and it will be used in this paper. 

Thus, a very promising and straightforward candidate for an improved dynamical 
approach of dislocations including scale effects is the translational gauge theory of dislo- 
cations [18, 19, 21, 24]. In the gauge theory of dislocations, dislocations arise naturally as 
a consequence of broken translational symmetry and their existence is not required to be 
postulated a priori. Moreover, such a theory uses the field theoretical framework which 
is well accepted in theoretical physics. From the point of view of theoretical physics, 
the gauge field approach of dislocations [18, 19, 21, 24] removes the singularity at the 
core of a moving dislocation at subsonic as well as at supersonic speed [25], the gauge 
theoretical approach in essence being reminiscent of nonlocality. Up to now, no solution 
for a nonuniformly moving dislocation has been given in the dislocation gauge theory or 
another generalized theory of elasticity. The aim of this paper is the examination of a 
nonuniformly moving screw dislocation, for the first time, in the framework of dislocation 
gauge theory. We will calculate the retarded expressions of the clastic fields as well as of 
the dislocation density and dislocation current tensors. This paper provides new insights 
on the change of the dislocation core structure while the dislocation is moving at sound 
wave speed. 

2 Gauge theory of dislocations 

In this section, wc bricfiy review the gauge theory of dislocations in the form given by Lazar 
and Anastassiadis [21], and Lazar [25]. 

In dislocation dynamics, the following state quantities of dislocations are of impor- 
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tance^ 



Tijk f^ik,j (^i'j,ki ^ij ~^ l^iji (1) 

which are called the dislocation density tensor and the dislocation current tensor, respec- 
tively. They are kinematical quantities of dislocations and they are given in terms of the 
incompatible elastic distortion tensor j3ij and incompatible physical velocity of the mate- 
rial continuum Vi. Their dimensions are: — 1, [vi\ — length/time, [Tj^-fe] = 1/length 
and [lij] — 1/time. The dislocation density and the dislocation current tensors fulfill the 
translational Bianchi identities 

^jkiTijk,i = 0, Tijk + hj,k — hk,j — 0. (2) 

In the dynamical translation gauge theory of dislocations, the Lagrangian density is 
of the bilinear form 

11 11 

C^T - W ^ -PiVi + - Dijiij - - aijPij - - HijkTijk. (3) 

The canonical conjugate quantities (response quantities) are defined by 

_dC n H -_9_^ 

^''-~d^i' W '''~di~j' ^&r~k' ^' 

where Pi, aij, Dij, and Hijk are the linear momentum vector, the force stress tensor, 
the dislocation momentum flux tensor, and the pseudomoment stress tensor^, respec- 

SI 

tively. They have the dimensions: [pj] = momentum/ (length)^ = Ns/m^, [D^j] — 

momentum/ (length)^ = Ns/m^, [aij] — force/ (length)^ = Pa and [Hijk] — force/length = 
N/m. 

The Euler-Lagrange equations derived from the total Lagrangian density C — C{vi, (3ij, lij, Tijk) 
are given by 

We add to £ a so-called null Lagrangian, jCn — CTijPij — P^Vi, with the 'background' 
stress cr^j and the 'background' momentum as external source fields, which satisfy: 
p^ — a^j j = 0. In terms of the canonical conjugate quantities (4), Eqs. (5) and (6) take 
the form 

Dij J + Pi— p^, (momentum balance of dislocations), (7) 
Dij + Hijk,k + CTij — (T°-, (stress balance of dislocations). (8) 



^We use the usual notations: (iij^k '■= dk(3ij and Pij := dtPij. 

^The moment stress tensor Tijk = —Tjik can be obtained from the pseudomoment stress tensor: 
Tijk = -H[ij]k (see [24, 26]). 
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Eqs. (7) and (8) represent the dynamical equations for the balance of dislocations. Eq. (7) 
represents the momentum balance law of dislocations, where the physical momentum is 
the source of the dislocation momentum flux. Eq. (8) represents the stress balance of 
dislocations. The force stress and the time derivative of the dislocation momentum flux 
are the sources of the pseudomoment stress. It can be seen in Eqs. (7) and (8) that the 
sources of the dislocation momentum flux and the pseudomoment stress tensors are the 
so-called effective momentum vector (pi—pf) and the effective force stress tensor (o"jj — a° ) 
(see, e.g., [19]). These effective flelds, which are the difference between the physical fleld 
and the background fleld, drive the dislocation flelds. 

The conservation of linear momentum appears as an integr ability condition from (7) 
and (8). It reads 

Pi — CTijj — 0, (force balance of elasticity), (9) 

where the time-derivative of the physical momentum vector is the source of the force 
stress tensor. 

The linear, isotropic constitutive relations for the momentum, the force stress, the 



dislocation momentum flux and the pseudomoment stress are 

Pi = PVi: (10) 

(7y = X5ijpkk + fJ-iPij + Pji) + 7(Ai - Pji), (11) 

Dij = diSijlkk + d2{Iij + Iji) + d2,{Iij - Iji), (12) 

Hijk — ciTijk + C2{Tjki — Tkji) + cz{5ijTiik — SikTuj), (13) 



where p is the mass density. Here iJi,X,^ are the elastic stiffness parameters, Ci, 02,03 
denote the resistivity parameters associated with dislocations and ^1,^2,1^3 are the iner- 
tia terms associated with dislocation currents. The nine material parameters have the 
dimensions: [fi,X,'-f] = force/ (length) ^ = Pa, [01,02,03] = force = N and [(ii, (^2, c^a] = 

SI 

mass/length = kg/m. 

The requirement of non-negativity of the energy (material stability) E — T + W > 
leads to the conditions of semi-positiveness of the constitutive moduli. Particularly, the 
constitutive moduh have to fulflU the following conditions [21] 

P>0, 

//>0, 7>0, 3A + 2//>0, 

d2 > 0, 4 > 0, 3di + 2d2 > 0, 

ci - 02 > 0, oi + 2o2 > 0, oi - 02 + 203 > 0. (14) 

If we substitute the constitutive equations (12) and (13) in the equations (7) and (8) 
and use the deflnitions (1), we flnd 

di0jj,i - Vjji) + (^2 + d3)0ijj - Vijj) + (^2 - d3)0jij - Vjji) +pi= (15) 
di6ij{f3kk - Vk,k) + + d3){(5ij - Vij) + (^2 - d3){/3ji - Vj^i) 

+ Ci{f3ikjk — (3ij,kk) + C2{l^jiMk — Pjk,ik + Pkj,ik — Pkijk) + C3 [Sij{f3ik,lk — l^ll,kk) 

+ {I3kk,ji - I3kj,ki)\ + (^i'j = ^ij' (16) 

which are a coupled system of partial differential equations for the flelds v and (3. 
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3 Equations of motion of a screw dislocation 



We now derive the equations of motion for an arbitrarily moving screw dislocation. We 
consider an infinitely long screw dislocation parallel to the z-axis and traveling in the 
xy-plane. The symmetry of such a straight screw dislocation leaves only the following 
non-vanishing components of the physical velocity vector and elastic distortion tensor 
(see, e.g., [3, 4]): Vz, /3zx, Pzy, and for the dislocation density and dislocation current 
tensors: Tzxy, hx, hy The equations of motion of a moving screw dislocation read 

(^2 + ds){/3zx,x + /3zy,y " Av^) + pVz 
zy,xy f3zx,yy zx 
{d2 - ds){/3^x - Vz,x) + C2{l3zx,yy " l3zy,xy) + (/^ " l)/3zx 
{d2 + d^){l3zy - Vz,y) + Ci{j3^x 

{d2 - d3)0zy - Vz,y) + C2{/3zy,xx " Pzx,xy) + (A« - ^)/3zy 

where A = d^x + dyy. In addition, the condition (9) reads now 

(/X + -f){/3zx,x + Pzy,y) = fy^z- 

Prom the equations (18)-(21) we obtain the following relations 

Cl _ C2 

d2 + dz _d2- ds 
and we may introduce the quantities 

/i + 7 
2 _d2 + d^ 

"-'-^^ 
9 d2 -\- ds 

= ■ 

yW + T 

Here li and Li are the 'static' and 'dynamic' characteristic length scales and Ti is the 
characteristic time scale of the anti-plane strain problem. Due to the conditions (14), 
they are non-negative: > 0, Li > 0, ri > 0. The velocity of elastic shear waves is 
defined in terms of the 'dynamic' length scale Li and the time scale Ti: 

4 = ^-^. (28) 

In a similar way, we introduce the following transversal gauge-theoretical velocity defined 
in terms of £i and Ti : 





(17) 


(/^ + 7)/3L, 


(18) 


- 7)/9L, 


(19) 




(20) 


- l)Pzy. 


(21) 








(23) 




(24) 




(25) 




(26) 




(27) 



and we find the relation 



2 „2 

(30) 



^1 _ Q-T 



Applying Eq. (22), the equations of motion (17)-(21) can be given in the form 

tIv, - LjAv, + v,^ (31) 

r2/3,, - £2 A/3,, - (l - ^) t;,,, + = /^o,, (32) 

r,^/3,, - A^,, - (^1 - 1) + = (3%. (33) 

If we assume that the following relation is vahd (see also [25]) 

Li = iu (34) 

then the field equations (31)-(33) decouple to Klein-Gordon equations. Equations (30) 
and (34) give the relation ax = ct, that means that we have only one characteristic 
velocity ct under this assumption. Particularly, the elastic fields fulfill the uncoupled 
Klein-Gordon equations 

[l + ilDrr]v, = vl (35) 

[l + £?nT]/3.. = /9L, (36) 
[l + £?nT]/3., = /3°„ (37) 

with the following (1 -|- 2)-dimensional d'Alembert operator (wave operator) 

n^^^du-A. (38) 

In addition, using (1) we may derive from (35)-(37) inhomogeneous Klein-Gordon equa- 
tions for the non- vanishing components T^xy, hx, hy 

[i + £?nT]7;., = T,%, (39) 
[i + £?nT]/.. = /L (40) 

[l+^ln^]hy = I%- (41) 

In field theory, Klein-Gordon equations describe massive fields (see, e.g., [27]). That 
means that a dislocation is a massive gauge field. 

Prom the condition (34), we find for the inertia term of a screw dislocation 

d2 + d^^^^pll (42) 
that it is given in terms of the characteristic length scale £i . 
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4 Nonuniformly moving screw dislocation 



Now we consider a nonuniformly moving screw dislocation at the position {C,{t),ri{t)) at 
time t. The dislocation line and the Burgers vector b — bz are parallel to the z-axis. The 
dislocation velocity has components: Vx = C{t), Vy = f)(t). 

At first, we want to find the gauge theoretical solutions of the dislocation density 
and the dislocation currents for a nonuniformly moving screw dislocation. Thus, we have 
to solve the equations (39)-(41), where the right-hand sides are given by the following 
sources 

T^xy = b5ix-mSiy-v{t)), (43) 

I% = -VxTl^. (45) 

Wc consider the situation where the source terms have acted after initial quiescence at 
t — > — oo. The solutions for the dislocation density and the dislocation currents are the 
convolution integrals^ 

(46) 
(47) 
(48) 

where the Green function of the (1 + 2)-dimensional Klein-Gordon equation is defined by 

[l + £?nT]G^^ = <5(t)(5(x-)5(y). (49) 
The Green function of the Klein-Gordon equation in (1-1-2) dimensions is given by [28] 

.KG _ 1 H{t- t/ct) _ / 74^2-^2'^ „2 _ ^2 , ..2 



T 

zxy 




zxy 


Izx 


^^KG 


zx^ 


Izy 




zy-i 



' ' cos -\, T'^x' + y\ (50) 

where H is the Heaviside step function. Substituting (43)-(45) and (50) into (46)-(48), 
the spatial integration can be performed to give the results 

1 /ct5't 



with the notations 



X^X-i{t'), y = y-r){t'), t = t-t', R'=X^+f, 

S'i^t^ tT^^t- — . (54) 



^The convolution is defined by: T^^^yix, t) = G^^ * T°^y = At' J dx'G™{x, t, x', t') T^^y{x' , t'). 
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Note that is called retarded time, which is the root of the equation = and is 
less than t, and (x, y) is the distance between the field point {x, y) and the position of 
the dislocation {i,r]). Here, ^ and rj are the positions in the x and y directions of the 
dislocation at time t', T4 and Vy are the velocity components of the dislocation at the 
same time t'. The retarded time tx is the time before t, when the dislocation caused an 
excitation of the dislocation field, which moves from to {x,y) in the time R/c'^. 

In addition, t — t-^ = R/ct is the time when the dislocation fields move from {^,,1]) to 
{x,y) with velocity ct- Thus, from each point x', the dislocation fields (51)-(53) draw 
contributions emitted at all times t' from — oo up to tT- That is the reason why Eshelby 
[5] said: 'The dislocation is haunted by its past'. 



4.1 Solutions for the elastic fields in terms of potential functions 

If we multiply Eqs. (35)-(37) with Dt and use the 'classical' result [6, 3], we obtain 

[i + £?nT]nT^. = /L,. + /l„ (55) 

[1 + ilDT] Dt^.. = T^.y,y + (56) 
[1 + ^^Dt] DT^y = r°^.,. + (57) 

as a set of fourth-order partial differential equations. Eqs. (55)-(57) have the two- 
dimensional form of the Bopp-Podolsky equations [29, 30] (see also [28]) in generalized 
electrodynamics, introduced by Bopp and Podolsky in order to avoid singularities in elec- 
trodynamics. As source terms only the classical dislocation density T^^y and dislocation 
currents 1°^, I^y are acting. 

Following Eshelby [6] , we express the elastic velocity and the elastic distortion in terms 
of 'potential' functions (F, A^, Ay): 

'^z — ^y,x ~ ^x,y^ (58) 
Pzx — F,y H 2" ^i" (59) 

Cr|-i 



^zy^-F,,-—A,. (60) 
Then, if A^, Ay and F are chosen to satisfy the subsidiary condition 

^.,. + ^,y + i^ = 0, (61) 

the field equations (55)-(57) become equivalent to the following inhomogeneous Bopp- 
Podolsky equations 

[i^iln^]n^F^Tly, (62) 

[i + i\n^]n^A,^-ily, (63) 

[i^iln^]n^Ay^il. (64) 
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The solution of these equations for an infinite medium is defined by 



± \J J. , 

= -G^^ * 1%, 



(65) 
(66) 
(67) 



where G^^ is the Green function of the Bopp-Podolsky equation (or wave-Klein-Gordon 
equation) 

[1 + ^^Dt] DtG^^ = 5{t)5{x)5{y). (68) 
As in generalized electrodynamics, we solve (68) with the help of two fields 

[1 + ^^Dt] G^'' = , DtG^ = 5{t)5{x)5{y) (69) 

and we get 

^BP^^n_^2^KG^ (70) 

where the first field is the Green function of the two-dimensional wave equation [7, 31] 



1 H(t-r/c-r) 



27r ^t^-ry4 



(71) 



and the second one is the Green function of the Klein-Gordon equation (50). Finally, we 
obtain for the Green function of the (1 -|- 2)-dimensional Bopp-Podolsky equation 



.BP_ 1 H{t-r/c^) 



27r ^t^-ryc'^ 



cos 



(72) 



Substituting from (43)-(45) for the case of a single dislocation, the spatial integrations 
can be performed to give the results for the potential functions of a screw dislocation 



b /"'t 1 
F = — / dt'^ 

dt 



— 

27r 



cos 



ctS^ 



b 

27r / 



, Vy{t') 



iT 

dt 



COS 



COS 



(73) 
(74) 
(75) 



In analogy to electromagnetic field theory, we can call (73)~(75) the (two-dimensional) 
'Lienard-Wiechert potentials' of a screw dislocation in the gauge theory of dislocations. 
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4.2 Radiation part of the elastic fields 



On the other hand, using (44) and (45), we want to solve directly the equations (55)-(57). 
Consequently, we have to solve the equations 



' zxy,x 



V T° 



2 zxy,xi 



zxy,yi 



(76) 
(77) 

(78) 



where the right-hand side is given by (43). It can be seen that only the elastic distortions 
contain acceleration parts. The elastic velocity does not depend on dislocation accelera- 
tion. Using the Green function (72), we obtain the result for the elastic distortions 



27rc| 



*T y (f) 



d_ 
dx 



oo 

iT 



1 — COS 



dt 



, Vx{t')Vy{t') 



f^zy — 



oo 

■oo 



dt 



St 

, Vxit') 



d_ 

dy 



tT 



dt 



, V,{t')Vy{t' 



1 — COS 



1 — COS 



d 
dy 



dt 



, 4 - y'{t' 



1 — cos 



+ 



dx 



dt 



1 — cos 



1 — cos 



T 



(79) 
(80) 



The solution of Eq. (76) for the elastic velocity has the same expression as given in Eq. (58) 
with (74) and (75). It can be seen that the fields given by (79) and (80) consist of two 
parts: 

(a) Field depending on the dislocation velocities 14 and Vy alone and proportional to 
spatial derivatives - velocity field 

(b) Field depending on the dislocation accelerations T4 and Vy and proportional to 
1/S't - acceleration field or radiation field. 

All our results depend on ix which is a characteristic length scale of the gauge theory 
of dislocations. 
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